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Abstract 

We study the variance and the Laplace transform of the probabihty law of linear 
eigenvalue statistics of unitary invariant Matrix Models of n x n Hermitian matrices as 
n oo. Assuming that the test function of statistics is smooth enough and using the 
asymptotic formulas by Deift et al for orthogonal polynomials with varying weights, 
we show first that if the support of the Density of States of the model consists of g > 2 
intervals, then in the global regime the variance of statistics is a quasiperiodic function 
of n as n — > oo generically in the potential, determining the model. We show next that 
the exponent of the Laplace transform of the probability law is not in general 1/2 x 
variance, as it should be if the Central Limit Theorem would be valid, and we find the 
asymptotic form of the Laplace transform of the probability law in certain cases. 
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1 Introduction 



Random matrix theory deals mostly with eigenvalue distributions of various ensembles of 
n X n random matrices as ri oo. Many questions of this branch of the theory can be 
formulated in terms of the Eigenvalue Counting Measure A/'n, A/'n(A„) is the number of 
eigenvalues in a given part A„ of the spectrum of a matrix M„. In a particular case of nx n 
Hermitian matrices eigenvalues are real and we can write for A„ C M: 

n 

A4(A„):=tl{A;"^GA„, / = l,...,n} = 5^XA„(A;")), (1.1) 

1=1 

where xa is the indicator of A. A more general object is a linear eigenvalue statistic, defined 

as 

n „ 

MniVn] := J^Vni^n = Tr (^„(M„) = / ^n{XWn{dX) (1.2) 
1=1 

for a certain test function 

V7„:M-*C (1.3) 

It is known that in many cases there exists a scaling of matrix entries (i.e. a choice of the 
scale of the spectral axis), such that for a sufficiently big class of n- independent intervals in 
fll.l|l (test functions in ()1.2|l ) the Normalized Counting Measure of eigenvalues 

Nn=Arn/n (1.4) 

converges weakly to a non-random measure A^, known as the Integrated Density of States 
measure (IDS) of the ensemble 

Nn N. (1.5) 

The corresponding scale (asymptotic regime) is called the global (or macroscopic). The 
convergence is either in probability or even with probability 1. We refer the reader to works 
[ini El IHO] , where this fact is proved and discussed for two most widely studied classes of 
random matrix ensembles, the Wigner Ensembles (independent or weakly dependent entries) 
and the Matrix Models (invariant matrix probability laws). Analogous facts are also known 
for many other ensembles. 

Besides, in many cases the measure A^ possesses a bounded and continuous density [ini 

El: 

N{d\) = p{\)dX, (1.6) 

We call p the Density of States (DOS) of ensemble. 

These results can be viewed as analogs of the Law of Large Numbers of probability 
theory. Hence, a natural next question, important also for applications, concerns the limiting 
probability law of fluctuation of the Normalized Counting Measures of eigenvalues or their 
linear statistics, i.e. an analog of the Central Limit Theorem of the probability theory. The 
question is not going to be completely trivial, because eigenvalues of random matrices are 
strongly dependent, and, as a result, the variance of the linear statistics ()1.2|) with a test 
function does not grow with n ( see e.g. formulas ()2.36|) - ()2.37|) below). Nevertheless, it was 
found that in a variety of cases fluctuations of various spectral characteristics of eigenvalues 
of random matrix ensembles are asymptotically Gaussian ( see e.g. El UHl UHl UHl UH 
l2Ill23ll21ll2niEniEniEZIEHlEi)-In particular, for the global scale and a test function 
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in (|1.2j) this requires, roughly speaking, the same order of magnitude of the entries for the 
Wigner Ensembles (see for exact conditions, similar to the Lindeberg condition of 

the probability theory), and for the Matrix Models one needs to assume that the support 
of the IDS ()1.5|1 is a connected interval of the spectral axis j^Hj: supp = [a,b]. The last 
result was obtained by using variational techniques, introduced in the random matrix theory 
in paper in order to prove 

Being applicable to Matrix Models of all three symmetry classes of the random matrix 
theory (real symmetric, Hermitian and quaternion real matrices), the variational techniques 
were efficient so far in the study of fluctuations of eigenvalue statistics only in the case where 
suppiV = [a,b]. In this paper we consider only the Matrix Models of Hermitian matrices, 
but for a general case of a multi-interval support of A^: 

q 

a := supp N = \J [ai, k], q > 1. (1.7) 

1=1 

In this case we can use recent powerful results by Deift et al [T7] on asymptotics of special 
class of orthogonal polynomials. We find that if (p is real analytic, then the traditional 
Central Limit Theorem is not always valid for the case q > 2. In particular, the variance 
and the probability law oscillate in n as n ^ oo, hence their limiting form depend on a 
sequence rij ^ oo. Moreover, the limiting probability laws are not always Gaussian. 

The paper is organized as follows. In Section 2 we study the variance of linear eigenvalue 
statistics in the global regime, i.e., for n-independent if in ()1.2|) . confining ourselves mostly 
to the case of test functions ip. We find that the variance is quasi-periodic in n in general 
and its frequency module is determined by the charges 

A = iV([a,+i,oo)), / = l,...,g-l, (1.8) 

determined by the IDS of ()1.5j) and its support ()1.7|1 . Hence, the variance has no limit as 
n oo for q > 2 and its asymptotic forms are indexed by points of the subset ff^^ c T"^"^, 
which is the closure of limit points of the vectors 

{{(3^n},...,{P,.,n})eV-\ (1.9) 

where {An}, / = 1, g — 1 are the fractional parts of (3in, I = l,...,g — 1. This phenomenon 
has been already found in certain cases [3 E] , but we give its general description. 

In Section 3 we study the Laplace transform of the probability law of linear eigenvalue 
statistics (|1.2p in the global regime, passing to the limit along a subsequence 

mn,{x)}, {/3g_in,(x)}) ^ x G M^-\ (1.10) 

and confining ourselves to real analytic test functions. We give first a general formula for 
the corresponding limit. Since the formula is rather complex, we consider several particular 
cases, where we show that the exponent of the limiting (in the sense fll.lO|) ) Laplace transform 
is not quadratic in ip, hence the limiting law is not Gaussian (see formulas (jSUl), fj3.13|) . and 
fj3.19|) - ()3.20|) ). This has to be compared with results of paper |T3], according to which the 
limits of variance and the probability law are the same for all sequences nj oo (i.e., exist), 
and the hmiting probabihty law is Gaussian. 
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The random matrix theory deals also with two more asymptotic regimes in addition to 
the global one. Namely, if, having fixed the global scale allowing us to prove ()1.5|) . we set in 

= ^ii\~ Ao)n") , < « < 1, Ao G supp N, (1.11) 
where if is n-independent, we obtain the intermediate regime, and if 

(^„(A) = (^((A - Ao)n), Ao G supp AT (1.12) 

then we have the local (or microscopic) regime. In Section 4 we discuss the form of variance 
and the validity of the CLT in these regimes. 

In Appendix A.l we find the asymptotic form of the variance of the Eigenvalue Counting 
Measure (jLip (corresponding to a piece- wise constant = xa in (EH)) for the Gaussian 
Unitary Ensemble (GUE). The variance proves to be O(logn) instead of being 0(1) in the 
case of a test function. This shows that the Lipschitz condition ()2.35|1 on test functions 
that we use to prove that the variance is bounded in n is pertinent (for a more precise 
condition see ^j). In Appendix A. 2 we compute the variational derivative of f3i in the case 
q = 2, which is used in Section 4, and discuss related topics. 

We note that a completely rigorous derivation of the results of this paper, especially those 
of Section 3, requires rather technical and tedious arguments. They will not be presented 
below. Rather, we confine ourselves to the presentation of results, their discussion, and 
outline of corresponding proofs. 

2 Variance of linear eigenvalue statistics 
2.1 Generalities 

Recall that unitary invariant Matrix Models are n x n Hermitian random matrices, defined 
by the probability law 

PnidMn) = Z-^ exp{-nTTV{Mn)}dMn, (2.1) 
where M„ = {Mjk}lk=„ = Wj, 

n 

dMn = Y[dMjj Y[ d'RMjkd'^Mjk, 

j=l i<j<k<n 

and : M — > is a measurable function, called the potential and such that 

\/ : M ^ E+, V{X) >{2 + 6) log |A|, |A| > L, (2.2) 

for some positive S and L. 

The hmit (|1.5p can be described as follows 123] • Consider the functional: 

Sy[m]=- / / log\X - fi\m{dX)m{dn) + / V{X)m{dX), (2.3) 
Jr Jr Jm 

where m is a non-negative unit measure. 

The variational problem, defined by (j2.3|) . goes back to Gauss and is called the minimum 
energy problem in the external field V. The unit measure minimizing (j2.3|) is called the 
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equilibrium measure in the external field V because of its evident electrostatic interpretation 
as the equilibrium distribution of linear charges on the ideal conductor occupying the axis 
R and confined by the external electric field of potential V. We stress that the respective 
minimizing procedure determines both the support a of the measure and the form of the 
measure. This should be compared with the variational problem of the theory of logarithmic 
potential, where the external field is absent but the support a is given (see fl2.5()|l ). The 
minimum energy problem in the external field ()2.3p arises in various domains of analysis and 
its applications (see recent book [31] for a rather complete account of results and references 
concerning the problem). 

The measure and its support a are uniquely determined by the Euler- Lagrange equa- 
tion of the variational problem p| IMj: 

K//(A) = F, Xea, (2.4) 

Veff{\) > F, \^a, (2.5) 

where 

V,ff{X) = ViX) - 2 /" log |A - fi\N{dfi), (2.6) 

J a 

and F is a constant (the Lagrange multiplier of the normalization condition A^(M) = 1). 

According to |9j (see also |2S1), if the potential V in 1)2.11) - ()2.2p satisfies the local Lipshitz 
condition 

\V{Xx) - y{\2)\ < C|Ai - Asr, |Ai|, IA2I < L, (2.7) 

valid for any L > and some positive C and 7, then (jl.Sp holds with probability 1, and 
is the minimizer of ()2.3|) . Moreover, if V satisfies the local Lipshitz condition, and the 
support ()1.7|) is a finite union of disjoint finite intervals, then ()1.6|) is valid [ini 123 and the 
Density of States can be written as 



p(A)=P(A)^i?,(A), AG a, (2.8) 
where -P(A) is a continuous function, 

^IrJx) = ^JbJz) , R,{z) = \{{z-ai){z-hi), (2.9) 

z=\+iQ i—^ 

and ^jRq{z) is the branch, determined by the condition: ^jRq{z) = z'^ + 0{z'^~^), z 00. 
To obtain these formulas, provided that the support ()1.7|) is given, we differentiate ()2.4|) . 
()2.6|) and obtain the singular integral equation 

v.p.|-— ^ = -^, AGcr. (2.10) 

Then the bounded solution of the equation has the form (|2.8p (see e.g. |28j) in which 

The endpoints of the support are rather complex functionals of the potential in general. 
Thus, it is of interest to mention a simple case |llj . 
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Let f : R — i> M be a monic polynomial of degree q with real coefficients. Assume that for 
some (7 > all zeros of f ^ — Ag are real and simple and set 

V(X) = ^M. (2.12) 

Then the DOS of the matrix model ()2.1|) with this potential is 

p{\)=''^\v\\)-Ag\''\\ea, (2.13) 
2'Kgq 

where 

a = {A G M : f^(A) < 4^}. (2.14) 
Besides, in this case we have for the charges ()1.8|) 

A = ^, / = l,...,g-l, (2.15) 
Q 

hence the set W'^ is {0, 1/g, (g - l)/g}. 

The case g = 1 corresponds to the Gaussian Unitary Ensemble and (j2.13p yields the 
semi-circle law by Wigner: 

In the case g = 2 and 

x;(A) = - m^, > 2^, (2.17) 

we have: 

a=[-h,-a][j[a,hl a=^m^- 2^, b=^m^ + 2^, (2.18) 

and 

... ,,,,, 

We will use in this paper the expressions for the variance of linear statistics (|1.2j) and for the 
Laplace transform of their probability law via special orthogonal polynomials. The technique 
dates back to works by Dyson, Gaudin, Mehta, and Wigner of the 60s (see e.g. |2I|). Namely, 
we have for the joint probability density of eigenvalues of ensemble ()2.1|) : 

Pn(Ai,...,A„) = (det{^]:\(A,)};_,^^)yn!, (2.20) 

where 

^g-ny/2pH^ (2.21) 

and 

{pI%>o (2.22) 
is the system of orthonormal polynomials with respect to the weight 

Wn = e-"^, (2.23) 
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so that 

[ e-"^Wp/")(A)Pir)(A)rfA = 5,,^, /,m = 0,l,... 
The polynomials satisfy the three - term recurrence relation for / = 0, 1, ...: 

r;"v;:i(A) + .;"v;"^(a) + rSwi(A) = a^-ha), r^-i = o, 

thereby determining a semi - infinite Jacobi matrix: 

Jj"k = rf^^j+i,k + 4"^5j,fc + j, A; = 0, 1, ... 

By using ()2.2()j) it can be shown that 





\v{\i)-v{^2)VKi{\A2)dXid\, 

2 

V„(Ai, A2)rfAi(iA2, 



Av9 



AA 



where 



n-l 



i^n(Ai,A2) = 5^^}")(AiM^"HA2) 



(2.24) 
(2.25) 
(2.26) 

(2.27) 
(2.28) 

(2.29) 



1=0 



is known as the reproducing kernel of the system (j2.2H) . 

Alp Lp{Xi) - Lp{\2) 



AA 



Ai — A2 



and 



V„(Ai,A2) = fraf^W(A0^a(A2) -^a(A0^i"HA2 



2. 



(2.30) 
(2.31) 



Note that in passing from ()2.27p to ()2.28|) we used the Christoffel-Darboux formula 

in) ^i"HA0^a(A2) - V^a(Ai)4"HA2) 
Ai — A2 



(2.32) 



It is of interest to have a spectral-theoretic interpretation of the above formulas. Notice that 
J*^"^ of ()2.26|) determines a selfadjoint operator, that we denote again J("). It acts in P{Z+), 
and the matrix 

S^-\dX) = e(")(A)rfA, eW(A) = {e^J{X)}rm=o, e!:^(A) = W"Ha)^L"^(A) (2.33) 
is its resolution of identity [T]. This allows us to write p.31|) in the form 

V.(Ai,A2) = (rt\)^(ea(A0ei"i,„_i(A2) (2.34) 
+ eS(A2)ei"2,,_i(Ai) -2e2i,„(A0ei':Li(A2))/2. 

Assume now that (f is of the class and 

sup|(/?'(A)| < cx). (2.35) 



for some positive C. It follows then from ^TM . (IT^ - (|T!TT|l that 

Var{MM} < {r^:li? fsup ■ (2.36) 

It can be shown that if the potential satisfies the local Lipshitz condition ()2.7j) . then the 
coefficients r^nli are bounded in n. We conclude that the variance of a linear statistic ()1.2j) 
is bounded in n if the test function satisfies ()2.35j) : 

Var{Ar„[v9]} < Const f sup |v9'(A)|) . (2.37) 

This has to be compared with the well known fact of probability theory, according to which 
the variance of a linear statistics of i.i.d. random variables is 0{n) for any test function. 

Notice that the condition ()2.35|) for ()2.37|) to be valid is pertinent. Indeed, it is shown in 
Appendix A.l that for the Gaussian Unitary Ensemble ()2.16|) and for the piece-wise constant 
V'n = Xa in ()1-2|1 . i.e., for the variance is O(logn), as n — >• oo (see formulas (jA.ljl . 

fpnn|i - dHH)). 

2.2 Asymptotic behavior of ?/'["^(A) 

We present now powerful asymptotic formulas for orthonormal functions of ()2.21|) - 
()2.24|1 due to Deift et al [Hj , valid in the case of a real analytic potential V in ()2.1|1 . We 
give a bit different their form that will be more convenient below. Our form is reminiscent 
of standard semi-classical formulas and is similar to that, given in [3 IE]- Notice that |Hj 
contains another, although heuristic, derivation of the asymptotic formulas as well as an 
interesting heuristic explanation of the quasiperiodicity of formulas for q > 2, based on ideas 
of statistical mechanics and quantum field theory. 

Assume that V is real analytic. Then the support p.7j) of is a union of g < oo finite 
disjoint intervals jT2j. The function P in ()2.8j) . ()2.11|1 is also real analytic. Following jT2j we 
say that V is regular if the inequality ()2.5p is strict and P is strictly positive on a. Hence 
the DOS p of ()1.6|) is strictly positive inside a and vanishes precisely like a square root at 
the endpoints. Denote 

Ar(A) = A^([A,oo)). (2.38) 

According to [Hj, in the regular case there exist functions dn{X), and 7n(A) such that if A 
belongs to the interior of the support p.7p . then 

V'i")(A) = (2rf„(A))i/2cos (7™Ar(A) +7„(A)) +0 (n-i) , n ^ oo. (2.39) 

Moreover, dn{X) and 7n(A) depend on n via the vector n/3, where (3 = {Pi, .., /?q_i) is given by 
(|1.8p . This means that there exist rz-independent continuous functions T> : a x T^^^ ]R_,_, 
and ^ : a X T*?-! ^ M such that 

d„(A) = P(A, nf3), 7(A) = ^^(A, nf3). (2.40) 

If A belongs to the exterior of cr, then il)^n\X) decays exponentially in n as — > oo. 

Similar asymptotic formulas are valid for coefficients of the Jacobi matrix J^"^ of ()2.26|) . 
i.e., there exist n- independent continuous functions TZ : T"^"^ M+ and S : T'^~^ M such 
that 

ri"_\ = 7^(r^/5) + O {n-^) , = S{nl3) + O {n-^) , n ^ oo. (2.41) 
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The functions V, Q, TZ, and S can be expressed via the Riemann theta-function, associated in 
the standard way with two-sheeted Riemann surface obtained by gluing together two copies 
of the complex plane slit along the gaps (61,02), (6g-i,ag), (6g,ai) of the support of the 
measure N, the last gap goes through the infinity. 

The components of the vector (3 = {(3i}'^~l are rationally independent generically in V, 
thus the functions V{X,n(3),Q{\,nP),7l{nP), and Sinjd) are quasiperiodic in n in general. 

We will also need asymptotic formulas for ^/'^'^^ as n — > cxd and k = 0{1) (in particular, 
we need the case /c = — 1 in ()2.31|) ). They can be extracted from [T7] (see also [8J)), but 
it will be convenient to present a different derivation, because the corresponding argument 
will be used in our analysis of limiting laws for linear eigenvalue statistics. To this end we 
replace a regular potential V in ()2.2()|1 by V/g, g > 0, introducing explicitly the amplitude 
of the potential. Then the quantities of asymptotic formulas fl2.38j) - ()2.41|1 will depend on 
g, and it follows from the results of [T71I2E] that these quantities will be continuous functions 
of 5^ in a certain neighborhood of 5^ = 1, provided that the support (jl.Tj) for (7 = 1 consists 
of q disjoint intervals. Consider now r^^^f^_i{g). Taking into account that the origin of the 
super- index n in the above formulas is the factor n in front of V in ()2.23|) . we can write 

V , V 



{n + k) (2.42) 
^(1 + k/n) 

In other words, to obtain r^n+k-iid) ^^"^ "^^n+ki^id) have to make the change 



g^g + -g (2.43) 

n 

in the inverse amplitude of the potential. We obtain in view of (j2.41|) for k = o{n): 

r^l-M = rit,-\((l + fc/n)^) (2.44) 
~ n{{l + k/n)g,{n + k)l3{{l + k/n)g)) 
- T^{9,n/3{g) + ka{g)) , 

where 

a{g) = igP{g))\ (2.45) 

and the symbol "~" denotes here and below the leading term(s) of the corresponding l.h.s. 
as — > cxD. 

We have an analogous formula for ()2.38|) : 

(n + k)N (A, (1 + k/n)g) ~ nN{X, g) + kv{\, g), 

where 

i^{X,g) = j-^{gN{X,g)). (2.46) 

By using these formulas, we can write for any fixed k (in fact k = o{n)): 

^l^^:l{\g) ^ {2V{X,g,n(3 + ka)y/' (2.47) 

iN{X, g) + 7rA;z/(A, g) + ^(A, g, nP + ka)^ . 



X cos 7m 
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Applying to ()2.44|) and its analog for s^^^lf^ the limiting procedure ()l.in|l . we obtain the 
coefficients 

rk-i{x) = n{x + ka), Sk{x) = S{x + ka), x G T''-\ k e Z. (2.48) 

They determine a family of the double infinite Jacobi matrices J{x), x G T''^^: 

{J{x)i))k = rk{x)i)k+i + Sk{x)i)k + rk-i{x)i)k-i, k e Z, (2.49) 

that can be viewed as a quasiperiodic operator, acting in 1^{Z). 
Consider now the functional 

S[m] = — log \X — fi\m{dX)m{d^), (2.50) 

defined on unit non- negative measures, whose support is contained in a. This is the standard 
variational problem of the potential theory i34|. Denote v the unique minimizer of ()2.50|) . 
Then, according to [12^ (see also [1111211), the non- increasing function u{\) = z/([A, oo)) 
(cf (I2.38|) ) coincides with the function, defined in (j2.46|) . Moreover, according to the 
measure z/ is the Integrated Density of States (IDS) measure of J{x) (see [32] for a definition 
of the IDS measure in a general setting of ergodic operators), the support ()1.7|) is spectrum 
of J(x), and (cf (Oil ) 

ai = u{[ai, oo)), / = 1, g - 1 (2.51) 

are the frequencies of quasi-periodic coefficients (j2.48|) of J{x). In other words, J{x) is a 
"finite band" Jacobi matrix, well known in spectral theory and integrable systems |41j . 

By using these facts and also taking into account that '?/'i"''(A) decays exponentially in 
n outside the support ^7], we can prove that for any continuous $ : M ^ C of a compact 
support we have in the limit (jl.lOp of ()2.33|) 

Ir / '^(^)<(.mn,(.)+"^(^)^^ = / meiM,^)dX, (2.52) 
where we have for A G cr: 

eim{X, x) = ifj^iX, x)iIj+{X, x) + V^;"(A, x)^-(A, x), (2.53) 



^+(A, x) = e'^'^^^^^UiX, x + la), U{X, x) = {V{X, x)/2)^/2gia(A,x)^ ^2.54) 

and ipf = ip^ . The above formula can also be written as 

^±(A, x) = ip^{X, X + T^a), (2.55) 

where T : T^^^ f''^^ is defined as Tx = x + a. This shows that \I/^(A, x) = {'ijjf{X, x)}i(zz 
is a generalized (quasi-Bloch) eigenfunction of J(x): 

J(x)^^(A,x) = A^^(A,x). 

In fact, {\E'^(A, x)}Aeo- form a complete system in /^(Z)). 

We refer the reader to [HI] for more details of this aspect of asymptotic formulas of [T7j . 
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2.3 Asymptotic behavior of variance 

We assume in this subsection that the test function in p.2|) is of the class and does not 
depend on n. Hence, function fl2.30|l is continuous in (Ai,A2). As a result, fast oscillating 
in n functions, entering ()2.47|1 . do not contribute to the limit, as was already in obtaining 
flTK^ - (12311). We have then from ^TMf . (jZHl, and (jTK^ - (jTK^ 

Var{Ar„[<^]} ~ V(n/5), (2.56) 

where 




AA 



V(Ai,A2,x)rfAirfA2, (2.57) 



V(Ai, A2, x) = 7l^{x) (eo,o(Ai, a;)e_i _i(A2, x) - e_i,o(Ai, a;)e_i,o(A2, x)) , (2.58) 
and eim{X,x) are given by (|2.53|) - (j2.54p . in particular 

eim{X,x)d\ = 5im. (2.59) 

Since the charges /?/, / = !,..., q—1 of (jl.Sp are continuous and non-constant functionals of the 
potential (see fZllin]), the leading term of variance is a quasi-periodic function generically 
in the potential. In particular, it has no limit as n — > 00. Its limiting points are indexed by 
the subset W^~^ C T*^"^, the closure of limiting points of the q — 1 dimensional vectors (jl.9p . 
W^^^ is T''^^ generically in V , but it can also be a proper subset of T^^^ (see, e.g. ()2.15|) ). 

The simplest case is g = 1 of a single interval support. Here H° is a point, there exist 
the limits [3 US] 

lim r^nlk-i = 1™ ■^i+fc = s, \/k e Z, (2.60) 

n— >oo n—>oo 

and the "limiting" Jacobi matrix J{x) of ()2.49|) has constant coefficients: 

Jjk = r5j+i^k + s5j^k + ^^i-i.fc- 

Placing the origin of the spectral axis at s, we obtain that A = 2rcos7rz/, ^/^^ are just the 
plane waves, 

a = [-2r,2r], 

and 

V{X) = = A G a. (2.61) 

This and general formula ()2.58|1 yield a version of (j236|l - (jT57jl in which the role of 
V(x, Ai, A2) plays 

This form of the variance was first found in physical papers |31 ^] and proved rigorously in 
[23- We see that in the single interval case the variance is universal, i.e., its functional form 
does not depend explicitly on the potential, the information on the potential being encoded 
in the unique parameter r of ()2.6()|1 . In particular, we have ()2.62|1 for the Gaussian Unitary 
Ensemble ^nkf [13]. 
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In the case (I2.12j) - (j2.15|) consists of q points (see (|2.15|) and the variance is a 

g-periodic function of n. Example ()2.17|) - ()2.19|) corresponds to the simplest non-trivial 
case q = 2, where Pi = 1/2, = {0, 1/2}, the matrix J(x) is 2-periodic, its coefficients are 



and the variance is asymptotically 2-periodic function in n [HI El ^ 

Var{A/;[v9]} ~ V(2)(n/2), 
where V^'^^x) is given by (jTKTjl in which V(Ai, A2, x), x eM^ = {0, 1/2} is 




(2.63) 



V(^)(Ai,A2,x) 



1 ^Xi^M 



(2.64) 



X 



((a^ - AiA2)(6' - A1A2) - (-l)2"a6(Ai - A2)^ ) , Ai, A2 G a, 



with 



i?2(A) = (A^ 



(2.65) 



and £a = 1 if A e (—6, —a) and ex = —1 if A G (a, 6). In fact, these formulas are valid for 
any real analytic and even potential, producing a symmetric two-interval support 



(see 13 E] ) ■ The general case of a two-interval and not necessarily symmetric support was 
analyzed in [H], where it was found that the variance can be expressed via the classical elliptic 
functions of Jacobi and Weierstrass. 

We conclude that a minimum modification of the limiting law of linear eigenvalue statis- 
tics in the case of a multi-interval support of the IDS, comparing with the case of i.i.d. 
random variables, could be a family of normal laws, indexed by the points of H^"^. We shall 
see below that this modification is not sufficient in certain cases. 

We remark in conclusion of this section that formulas ()2.56|) - ()2.58p allow us to charac- 
terize the universality classes of ensembles (|2.ip with respect to the variance in the global 
regime, i.e., the sets of ensembles (potentials), leading to the same asymptotic form of the 
variance of linear statistics in the regime. Namely, since the potential is present in ()2.56|) - 
()2.58|] only via the endpoints (ai,...,6g) of support and via the charges (/3i, /3g_i) of all 
but one intervals of the support, these parameters determine a universality class. Notice 
that the parameters are not necessarily independent. 

3 Limiting laws 

3.1 Laplace transform of the probability law of linear eigenvalue 
statistics 

In this subsection we obtain an expression for the Laplace transform of the probability law 
of linear eigenvalue statistics ()1.2|) via orthogonal polynomials. We consider here real- valued 
test functions : R ^ M. The Laplace transform is evidently 




(2.66) 




(3.1) 
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where Ey{...} denotes the expectation with respect to (|2.1|) (or (|2.2(jp ). determined by a 
given potential V, and 

Afn[v]=KM-'Ev{MM}- 

It is convenient to introduce the parameter s G [0, 1] and to consider the function 

= iogZ„M, s e [0,1]. 

It is easy to see that 



F„(0) = 0, F^(0) = -Ev^ l^nb]} 



and 

= Ey.,^ {^^M} - E^+^ {KM} ■= Varv^+^ 
This yields the following expression for the logarithm of (|3.1|) : 

logZ„[yp] = Fn{l) := F4^] = / (1 - s)VsLrv+^ {Afn^jds. (3.2) 

We mention that there exist another expression for the Laplace transform p.ip . It dates 
back to the Heine formulas in the theory of orthogonal polynomials (see e.g. ^Hj, Theorem 
2.1.1) and can be easily obtained from the Gram theorem: 



n 



j,k=l 



Z4^] = det / e-'^(^Vf^(A)V^f)(A)rfA 
= e^{-^"[^]>det(l-ir„,^), 

where 

(^(A) = ^(A)-E^{Ar„[<^]}, Ev{ArM} = n [ <^(A)Ey{iV„(dA)}, 
and Kn,ip is the integral operator, defined as 

(ir„,^/)(A)= / ir„(A,/i)(l-e-'^('^))/(/i)rf/i, AgR. 
Jr 

These formulas and their analogs for unitary matrices were used to prove various versions of 
the Central Limit Theorem (see e.g. [HI |2S1 121 EH EHl EHl HI)- 



3.2 Asymptotic behavior of the Laplace transform 

We will assume in this subsection that (f is real analytic. According to ()3.2|) . ()2.28|) . and 
fl2.3H) . we have to find the asymptotic form of iplC^ and V'i-i for the potential V + s(p/n. 
We have already seen in the previous section that adding terms of the order 0{n~^) to 
the potential we obtain non-trivial contributions to the asymptotic formulas because of fast 
oscillating in n functions in the r.h.s. of (j2.39p - (j2.41|) . etc. The 0{n'^) terms appeared 
there because of the passage n ^ n + k, leading to ()2.45p - ()2.48|) . In this case the terms are 
proportional to the potential, since we just change its amplitude: V — >• V{l — k/n) (see ()2.42p 
- ()2.43|) ). This required derivatives ()2.45|) and ()2.46|) of "frequencies" Pi, I = 1, ...,q — l, and 
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iV(A) of fast oscillating functions in (j2.39|) - (j2.4H) with respect to the inverse amplitude g 
of the potential. 

On the other hand, to find the asymptotic behavior of the Laplace transform, we have 
to add to the potential the term sip/n (see (13. 2|) ). Since ip V in general, this requires 
variational derivatives of frequencies with respect to potential, i.e., we have to add the term 
ecp to the potential, and find the derivative of Pi, I = 1, ...,q — 1, and A^(A) with respect to 
e at e = 0. 

Consider first the case q = 1, where the support of the IDS is a single interval. Here the 
dependence on x of functions V, G, 7^ and S of - pTHl is absent (see - JZH). 

Hence the term s(p/n is negligible in the limit n —>■ oo, because there are no fast oscillating 
in n functions in the asymptotics of ip^^+k' ^ ~ O5 "'"n^i obtain from 

(K^ and dsn: 

lim F^[if] = lim Var{Ar„[v5]}/2. (3.3) 

n^oo n— >oo 

Notice also that we have here the "genuine" limit as n — > 00, but not a sublimit (|l.l(Jj) along 
a subsequence. We conclude that the Central Limit Theorem is valid in this case. This was 
proved in |23| by the variational method and for a rather broad class of potentials and test 
functions (not necessarily real analytic). 

As is was shown in the previous section, the variance of a linear statistics with a test 
function has no limit as n ^ 00 if g > 2. Its sublimits are indexed by points of the "hull" 
jjg-i f9-i^ Hence we cannot expect the traditional CLT ()3.3j] . as in the case of g = 1. 
Rather this should to be a collection of the CLT, indexed by ff^^: 

lim F„^^(.)[<^]= lim Var{Ar„^(,)[v9]}/2 = V(x)/2, x G e''-\ (3.4) 

where {nj{x)} and V{x) are defined in ()1.10|) and ()2.57|) - ()2.58|) . We will call this the 
generalized CLT. 

We will show now that the generalized CLT is not always the case for q > 2. Recall that 
A^(A) and Pi are functionals of V and denote 



, / = l,...,g-l, N[^] = |-iV(A) 
£=0 '^^ 



(3.5) 



e=0 



the variational derivatives of Pi and A^(A) with respect to V. Pi[(p] and N[(p] are linear 
functionals of if and nonlinear functionals of V. It follows from that they are well defined 
if V is real analytic and regular and ip is real analytic and such that maxAgR |^(A) / ^^(A) | < 00. 

Arguing as in Section 2.2, we obtain that in this case is given by ()2.47p - ()2.48p 

with the replacement 

kai — >• kai + sPi[ip], nku — > Trku + nsN[ip]. 
Now, assuming (|1.1U|) and taking into account (j3.2j) . (j2.28|) . and (|2.34|) . we obtain from (|3.2|) : 

F[^]:= lim logZn^(,)[y^]= [ {l-s)V{x + sP[(p])ds, (3.6) 

where V is given by ()2.57j) . According to ()2.57j) V is a quadratic functional of (p. Hence the 
functional F[ip] is not quadratic in general, because of the presence of the term sP[(p] in the 
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argument of the integrand in (j3.6|) . In other words we have here a hmiting law in the sense 
of p.lU|) . but the law is not necessarily Gaussian. 

It seems that a general classification of possible cases is rather complex. We thus will give 
several examples, showing different cases of asymptotic behavior of the Laplace transform 
of the probability law of linear eigenvalue statistics. 

Consider first the case, where the test function is a multiple of the potential: 

ip{X) = tV{X), t G R. (3.7) 
Then (ITIHl) and the relation pilV] = - l3'{g)\g^^ yield: 

= PiitV] = -t - aiil)) , / = 1, g - 1, 

where a;, / = 1, g — 1 are defined in ()2.51|1 . Hence, if 

A(i) = aKi), / = i,...,g-i, (3.8) 

then the integrand in ()3.(j|l does not depend on s, and we obtain in view of (j2.5(i|l - ()2.57j] 
the generalized Central Limit Theorem ()3.4|1 . 

The equality ()3.8|) is valid for any potential of the form ()2.12|) with g = 1 and f ^ — 4 
having only simple and real zeros, because according to ^] 

A(1) = «K1) = ^, / = l,...,g-L 

q 

It is also valid for any even potential, having two equal local minima and one local maximum, 
which is high enough to produce a two-interval support (|2.66p . In this case (13. 8p results from 
the symmetry, implying 

f3^ = ai = 1/2 (3.9) 

(recall that in this case the vectors j3 of ()1.8|) and a of ()2.45p are one- dimensional: (3i = 
N{a), ai = z^(a)). 

In all these cases the limiting Jacobi matrix J of ()2.49|) is g-periodic (g = 2 in the case 
of (ESSl), see (ESni))- 

It can also be shown that we have the generalized Central Limit Theorem for potentials 
fl2.12p and (p = tv (here the limiting matrix J is also g-periodic). 

To demonstrate a possibility to have a non-Gaussian limiting law, we consider a simplest 
non-trivial case of even potential with the two-interval support ()2.66|) and of test function 

^(A) =tX, te M, (3.10) 

i.e., the case of "linear" linear statistic 

n 

t^Aj"^ =tTr M„. (3.11) 

1=1 

Since in this case Aip/AX of (j2.3Up is equal to t, it follows from (j2.58|) and (j2.59|) that 

V(x) = ^27^2(x), (3.12) 
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and then (|3.6|) implies that in the case (|3.1U|) (and for any support) we have for the exponent 
of the hmiting Laplace transform: 



ip{\)=t\ nj(x)^oo 



lim logZ„^,(,)[^] =e ! {I~s)n^x + s'f5y\ 



ds. (3.13) 



According to |8j|), it is possible to express the coefficient TZ^x), corresponding to the two- 
interval support, via the Jacobi elliptic function: 



n'{x) = + ^ cn\x + 1/2), (3.14) 

where cn(a;) = cn{2K{k)x\k), k"^ = Aab/{a + 6)^, K{k) is the elliptic integral of the first 
kind. In view of ()3.9|) the coefficient of ()2.48p is given by ()2.63|) : 

^fk\ b-i-D'^a 



ip{X)=tX 



and is 2-periodic (see also j2]). In view of (j3.12|) this implies that the variance of (j3.11|) is 
asymptotically 2-periodic in n: 

(^) = ^ - (-1)4 (3^15) 
Furthermore, it is shown in Appendix that 

Hence we obtain from ()3.14|) 

FM^ix)=tx = / - ^)^'(^ + (3-17) 

^0 

It follows from (j3.16p that uj is irrational generically in a, and b, hence TZ'^{x + su) is quasi 
periodic in s in these cases. Since TZ^ is 1-periodic and real analytic, we can write its Fourier 
series 

7e2(x) = ^c^e^™ (3.18) 
with fast decaying coefficients. Plugging (j3.16|) . and (|3.18|) in (|3.13|) . we obtain 

FM^ix)=tx = ^ - ^^'(^) - ^(^) + ^(^ + ^^)' (3-19) 

where 

A{x)= V TT^^^e^— . (3.20) 

^ ^ ^ (27r«mcj)2 ^ ^ 

We see that the logarithm of the limiting Laplace transform of the probability law of statistics 
(j3.11|) contains not only a multiple of t^/2, that would correspond to the CLT, but also a 
linear in t term, a constant in t term, and either quasi periodic (generically in a, 6, when lo 
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is irrational) or periodic (in special cases, where u is rational) function of t. Besides, while 
the variance of statistics ()3.1H) is ()3.15p in the limit p.lO|) . the coefficient in front of is 



Co = / TZ^{x)dx 
Jt 



b^ + a^ 



hence is not the variance ()3.15p . 

Notice also that in the case g > 3 we would have in ()3.20p the sum over Z''~^\{0}, and 
the expression 

p ■m:= p^rrii + ...pg_^mg_i, 

in the denominator, that can be arbitrary small for certain collections of (mi, mg_i) G 

Z^^^ and (3 := {(3i, l3g_i) G in the case, where the components of (3 are irrational. 
Hence, to make the series in ()3.20|) convergent, we have to assume that the components of 

(3 are sufficiently bad approximated by rationals (e.g. a Diophantine condition). 

According to Appendix, in a general case of a real analytic if and a two-interval support 

cr = [ai, bi] [J [a2, 62], -00 < ai < 61 < 02 < 62 < 00, (3.21) 

we have: 

where ^/RqJjj) is defined in fl2.9|l . 

/ = r = 2 

Jb, V(&2-/i)(a2-/i)(/i-&i)(/i-ai) {{b2-h)ia2-a^)y/^ 

and 

^2 _ (Q2 - bi){b2 - ai) 
ih - 6i)(a2 - oi) 
(see [20]; formula (3.149.4)). In the symmetric case ()2.66|) we have 



(3.23) 



2 .A- 2^ U^/f^" 



b + a \b + al a \a 

where the second equality results from the formula (1 + k)~^ K {2\/k / {1 + k)) = K{k) 
formula (8.126.3). It follows then that in this case and for ip = t\ ()3.22|) coincides with 

dSISI)- 

Combining (|3.6|) for g = 2 and (|3.22j) . we obtain the limiting form of the Laplace transform 
of the probability law of linear eigenvalue statistics for real analytic test functions. In fact, 
this form is also valid for larger classes of test functions, in particular, for bounded 
functions with bounded derivative. Indeed, for any such function ip there exists a sequence 
{ifk} of real analytic test functions, converging to ip in the metrics. Besides: 

(i) V of ()2.57j) is continuous in ip in the metrics and in x G T^; 

(ii) (3i of ()3.22j) is continuous in ip in the C metrics; 
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(iii) we have for the functional F„ of (|3.2|) 

\Fn[^i] - Fn[ip2]\ < C fsup - i^'S] , (3.24) 

Vagk / 

where C depend only on the potential and is finite if ()2.2p holds. To prove this inequality, 
implying the uniform in n the continuity of Fn in ip, we use an argument, similar to that 
proving ()H.2|1 . to obtain 

F[ifi] - F[if2] = / (1 - s)Varv+^,/nWn[fi - f2]}ds, 
Jo 

where v?, = sipi + (1 - s)ip2, s G [0, 1]. This and (jOTj) yield ^HM . 
Now we can write: 

for Fn of ()H.2|1 and F^ of (jH.fij) with explicitly indicated dependence on x G T^. Making the 
limit nj{x) —>■ oo and using ()H.24|1 . we obtain in view of the above result on real analytic 
test functions: 

lim sup \Fn^^^)[ip] - F^[ip]\ < C I sup - v^'fcl ) + l-F^M - F^[^k\\. 

Now the limit ipk —>■ ^ and above assertions (i) and (ii) yield the validity of our re- 
sults fl3.6|) for g = 2 and ()3.22|) for bounded test functions with bounded derivative, 
hence a possibility to have non-Gaussian limiting laws (the non-validity of the Central Limit 
Theorem) for linear eigenvalue statistics of ensembles 1)2.11) with these test functions. 

Formula ()3.22j) allows us to characterize the class of potentials and test functions for 
which the (generalized) Central Limit Theorem ()3.4|) is valid in the case of a general two- 
interval support ()3.21|) . Indeed, for any pair (V^, v?) for which the r.h.s. of ()3.22|) is zero 
the integrand of ()3.6|) does not depend on s and we have the generalized CLT ()3.4|) . In 

particular, in the symmetric case ()2.66p it follows from (|3.22|) that (3i [ip\ is zero if and only if 

^(/.)-^(-A.) ^^^^^ 



In particular, for an even potential of support ()2.66p and an even test function the gener- 
alized CLT is valid. 

One can view (p as an analog of external field in statistical mechanics. Hence, we can say 
that in this case an even external field " does not break the symmetry" . On the other hand 
a "generic" ip or an odd such that 



L 



is a "breaking symmetry field" and leads to a non-Gaussian limiting law. Its simplest case 
ip{X) = tX (pmUl is given by (nTTTl . 
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4 Intermediate and local regimes 

In this section we consider limiting laws of linear eigenvalue statistics for the test functions, 
given by (jl.llj) and (jl.l2j) with a function (p and Aq belonging to the interior of a. 

We begin again by calculating the asymptotic form of the variance in these cases. Chang- 
ing variables to 

Ai,2 = Ao + ti,2/ri°, < a < 1, (4.1) 

we obtain from ()2.28|) : 

Var{Ar4<^]}= / / ('^Vv„(Ao + ti/n",Ao + t2/n")rftirft2. (4.2) 

To find the asymptotic form of the r.h.s. we will use again ()2.H4|1 and ()2.HH|1 in which 
'0i+fc(A), = 0. — 1 is given by (j2.47j) with A = Aq + t/ra". Taking into account that V, N, 
and Q are smooth functions of A in a sufficiently small neighborhood of Aq, we can write 

V'it(A) ^ {2V{Xo,nf3 + ka)f' 

X cos (7™A^(Ao) + 7iku{Xo) - 7rp(Ao)n^-" + ^(Ao, np + kv)) , 

where p(A) = — A^'(A), and we do not indicate the dependence on g (in fact it suffices to 
consider the case g = 1). Plugging this into ()2.HH|1 and ()2.H4|1 . and omitting in the resulting 
integrand of the r.h.s. of ()4.2|) the fast oscillating terms, we obtain 



Var{Ar„[v?]} ^ fi(Ao,n/3) (4.3) 

rp(Ao)(ti - 




[^[ti) - (p{t2)) — dtidt2, 



where 



B{\,x) = 27i^n\x)V{X,x)V{X,x-a) (4.4) 
X sin^ (vrz/(A) + ^(A, x) - g{X, x - a)) . 

This leads to the following result in the local regime a = 1 and for the limit ()1.10|) : 

lim V8ir{Ma^(x)[<^n,{x)]} = B{Xo,x) (4.5) 

rij {x)—*oo 



It is known that in the local regime the variance of linear eigenvalue statistics has a universal 
limiting form J?! EH] 

lim Var{iV„K]} = / / (^(t,) - y^(t2))^ '^"' '^^^ ^t,cit2, (4.6) 



27r2(ti-t2)2 

in which all the information on the potential is encoded in p{Xq). This and ()4.5|) imply the 
identity 

B{X,x) = l, (4.7) 



19 



valid for all x G T^^^ and A, belonging to the interior of the a (a direct proof of the identity 
can be extracted from the proof of Lemma 6.1 of |[I7j ) . 

In the intermediate regime < a < 1 we still have a fast oscillating factor 



sin^ (ni-"7rp(Ao)(ti-t2)) 



in the integrand of ()4.3|) . Replacing the factor by its average 1/2, and using ()4.7p we obtain 
in this regime 

hm Var{iV„[^„]} = [ [ ^^^^^y^^^^dt.dt,. (4.8) 

As in the case of the local regime the limit here is the same for any subsequence fll.l()|l . hence 
we do not need to assume ()1.1()|1 . We conclude that the variance of linear statistics has a 
well defined limit in the intermediate regime as well. Moreover, ()4.8p is the "smoothed" 
version of variance ()4.6p in the local regime, since ()4.8|) is ()4.6|) in which the "oscillating" 
factor sin^(7rp(Ao)(ti — ^2)) is replaced by its average 1/2. 
Passing to the Fourier transform 

we can rewrite the r.h.s. of (14.81) in the form 

\k\(p{k)^{-k)dk 



I 



appearing in the continuous analog of the strong Szego theorem (see jHllIHlEZl and references 
therein) . 

The universality property of unitary invariant Matrix Models fTli IH^ implies that the 
Laplace transform of the probability law of linear eigenvalue statistics has the following 
limiting form in the local regime 

lim Zn[^„] = e2-^(^«)^(°) det(l - S^), (4.9) 



where S^, is the integral operator, defined as 



!5^i«^(l-e-<"))/M*, (4.10) 

7rp(Ao)(t -m) 

and we assume that in ()1.12j) is continuous and integrable on M. It is obvious that the 
logarithm of the r.h.s. of ()4.9|) is not quadratic in cp, hence the CLT is not valid in the local 
regime (see e.g. for related results). 

If, however, we take in the above formulas 

<f{t) = mt-to)6), 

where $ does not depend on 6 —>■ 0, and to G M, i.e., we assume that the test function in 
()4.9|) - ()4.1U|) is "slow varying", then it can be shown (see e.g. 0113 EH]) that the limit of 
the r.h.s. of as 5 ^ is the r.h.s. of divided by 2. 
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On the other hand, take as a test function in (|3.6|) 

ip{X) = mX-Xo)/6), (4.11) 

where $ does not depend on 5 — 0, and Aq belongs to the interior of the support of A^, i.e., 
assume that (p is "fast varying". Since the variational derivatives (linear functionals of ip) 

I = 1, g — 1 of ()3.5p can be written as 



= / kiXMX)dX, 

we have 



Jr 



Hence the term s(3 in the argument of the integrand of ()3.(jj) vanishes in the limit 5^0 and 
we obtain from ()2.57|) . changing variables to Ai,2 = Aq + 6ti^2- 

hm F[^] = iv(x), (4.12) 

5^oo <^(A)=${(A-Ao)/5) 2 

where 

V(x)^V(Ao,Ao,x)/ f ^l^'f^!f dt.dt.. 



Now it can be shown, by using (jZ^ . (jT^ . and (jO|) . that the r.h.s. of the last 

formula coincides with (j4.8j) . Hence the limit (j4.12j) coincides again with the r.h.s. of (j4.8p . 
divided by 2. 

The above suggests that the CLT is valid in the intermediate regime. This was indeed 
proved in several particular cases (see [HZ] and references therein). 



A Appendices 

A.l Variance of the Eigenvalue Counting Measure of the GUE 

It has been shown in Sections 2.1 and 2.3 that the variance of linear eigenvalue statistics is 
of the order 0(1) if the corresponding test function is of the class . Here we will argue 
that this condition is close to be optimal. To this end we will show that linear eigenvalue 
statistics for the GUE (|2.16p are of the order O(logn) if is the indicator xa of an interval 
A = (a, 6) G (— 2w,2w). For this choice of (/? the corresponding linear statistic is the 
Eigenvalue Counting Measure ()1.1|) . We will derive the asymptotic formula 

Var{Ar„(A)} = ^ logn (1 + o(l)), n ^ oo. (A.l) 

We have from (fTTfjl 

YariKiA)} = f dX [ K^iX, fi)dfi, (A.2) 
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where Kn is given by ()2.29|) . According to ()2.32|) we need asymptotics of ip^^'' for I = n — l,n 
to find the asymptotic behavior of ()A.2|1 . These follow from ()2.47j) (with the remainder o(l)), 
()2.6H) . and ()2.1fi|l or can be obtained from the well known Plancherel-Rotah formulas for 
the Hermite polynomials (see jlO], formula (8.22.12)): 

^it (^) = (2vry^Lg)V2 ^^«^^) + + ^^^)) ^1 + ^^1))' ^^-3) 

where A = 2^cos^, a{e) = 6- sin 2^/2, 0(6) = a{e)/2 - n/A. 
Besides, in this case jlOj 

= V^. (A.4) 

Var{Ar„(A)} = I, + !, + !, + I,, (A.5) 



We write flA.2jl as 
where 



/•oo /"b — a 

h = / c/x / Kl{b-y,b-y + x)dy, (A.6) 

pb—a px 

h = dx Kl{b-y,b-y + x)dy, 

Jo Jo 

and I3 and /4 can be obtained from Ii and I2 by replacing -ft'„(6 — y,b — y + x) by -ft'„(a + 
y,a + y-x). 

By using the ChristofFel-Darboux formula ()2.32|) we can write the integrands in Ji 2 as 

i&r^lix,y)x-^ 

where 

^nix,y) = fe")(A)^a(/.) - V^a(A)V^i")(/i)) _ ^ . (A.7) 
In view of ()A.4|) the integral Ji is bounded from above by 

gib -a)'' [ f^(")(A)V^a(^)-^a(A)^(")(/x))'rfArf/x = 2^/(6 -a)^ 



where in writing the last equality we used the orthonormality of the system {ipl""^} 

To find the asymptotic behavior of I2 we will use formula ()A.3|1 . However its direct use 
is impossible since the remainder term leads to the divergent integral in x. Thus, we write 
this integral as the sum of the two, over the intervals (0, A/n) and {A/n, b — a), where A is 
fixed. In the first integral we write ()A.7j) as 



«>0- 



X=b—y, pL=b—y+x 

In view of (jA.3|l i^^nlki ^ = 0, — 1 are of order 0(1). Besides, we have the relation [10] 

d ,{n) ^ \ iin) , ff^l\ , (n) 
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implying that -^i/j^^^. is 0{n), k = 0,-1. This yields the bound |\E'„(x,y)| < const ■ nx, 
according to which the first integral is 

"A/n 

const ■ / xdx = 0(1), n oo. 







In the integral over {A/n,b — a) we use the asymptotic formula ()2.47|) . neglecting, as it was 
done in the proof of the previous sections, the fast oscillating terms, and noticing that the 
remainder term o(l) in ()2.47|) yields the error o(logra). This and ()A.4|) for / = n lead to the 
asymptotic expression fcf I2.62p 



dx r , 4^ - A/i 

-T / dy 



+ 0(1) 

X=b-y, fj.=b-y+x 



The leading contribution to the integral is due the integral in x over an interval {A/n,e), 
where e > is small enough and n-independent, because the integral over (e, 5 — a) is n- 
independent, hence is 0(1), n oo. Then the condition \ = b — y, fi = b — y + x can be 
replaced by X = b — y, fj, = b — y. As a result the integral in y is asymptotically x, and we 
obtain ^ 

Ii + l2 = 7-^logn + 0(1), oo. 

The same contribution is due the sum J3 + J4 in ()A.5|) . hence we obtain ()A.1|) . 

Analogous formula is also known for unitary matrices. For this and corresponding Central 
Limit Theorem see EH US . 



Consider the covariance Cov{A/'„(Ai), A/'„(A2)} of the Eigenvalue Counting Measures 
(dUI) for Ai = (ai,6i) and A2 = (02,62), Ai,2 C [-2^,2^). It follows from (IT^ that 

(cf (1123)) 

Cov{Af4^,],Afn[^2]} := E{Af4ipi]Af4ip2]} - E{Ar„[<^i]}E{Ar„H} 

= 11 / U(Ai)-y.i(A2))(^2(Ai)-<^2(A2))i^^(Ai,A2)f/AirfA2. 
JrJr ^ ^ ^ ^ 

This yields for ipi^2 = XAi,2 

Cov{A/;(Ai),Ar„(A2)} = / [x(AinA2)xK(A,/i) - XAixA2(A,Ai)] Kl{X, fi)dXdi^. (A.8) 

The argument, proving ()A.1|) . allows us to find the leading terms of the covariance as n ^ 00 
in various cases. We will assume below for the sake of definiteness that ai < 61, 02 < 62, bi < 
62, and 61 — Oi < 62 — 02, and consider the following cases of the asymptotic behavior of 
flA.8|) (see for analogous results for unitary matrices), 
(i). Disjoint intervals (bi < 02): 

^ ^ dA / ,. d,, (A.9) 



4"-" J A, Ja, 

i.e., in this case the covariance is bounded as it was for the case of linear statistics, generated 
by functions in (ITirTI) and (E3T)D . 
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(ii). Touching (outside) intervals (02 = &i): 



27r 

(iii) Touching (inside) intervals (ai = a2, but bi < 62 

1 



2 logn. (A.IO) 



(iv) Embedded intervals (02 < cti, &i < ^2): 



logn. (A. 11) 



' ^ dA / rf/z. (A.12) 



(v) Intersecting intervals (02 < &i but ai < 02 and 61 < 62). In this case we can write the 
equality 

X(AinA2)xM - XA1XA2 = X(ai,a2)x{R\{ai,fe2)) X{ai,a2)x{bi,b2)^ 

(A.13) 

showing that the domain of integration in ()A.8|) does not contain the "dangerous" diagonal 
X = ^. Hence, the leading term of the covariance in this case is of the order 0(1), and the 
respective coefficient is given by the integral of the product of the integrand of ()A.12j) and 

of (jUHl). 

By using similar reasoning, it is possible to find the asymptotic form of the variance and 
the covariance of the Eigenvalue Counting Measure in other cases. Note that in the "regular" 
cases, where the covariance is bounded, (see ()2.56|) and ()A.9|) ). its leading term is additive in 
Ai and A2 (or in ipi and ip2), while in the "singular" cases, where Cov{A/'„(Ai), A/'„(A2)} = 
O(logn), its leading term is independent of Ai and A2. 



A. 2 Variational derivative of frequency in the two-interval case 



We derive here formula ()3.16|) . We will use the variational approach, based on the functional 

Write the minimum condition ()2.4j) . ()2.fi|l for V + e(p, and compute its derivative aX e = 0. 
This yields 

f{X) — 2 / log |A — fi\p{fi)d^ = const, A G a, (A. 14) 

J a 

where 

d 

and p is the density of the measure of Notice that the differentiation of the limits 

of integration in ()2.(ij) does not contribute to ()A.14j) . because p vanishes at each endpoints 
of the support according to ()2.8|) . 

The derivative of ()A.14|) in A is the singular integral equation (cf ()2.10|) ): 

r p_{pw ^ _m ^ ^ ^ 

y . /i - A 2 ' 
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The general solution of the equation in the case (|2.66p is [2B1 

CiA + C ^ 1 f ip'ifi)X2{fi)dfx 



X2iX) 27r2X2(A) "7, /i-A 

where 



X2iX) = -V^2(A), Vi?2(A) = ^/R^] 

z=X+iO 

Rliz) = {z^ — a^){z^ — IP') (see (12. 9p for q = 2) and R2{z) is the branch of the square root, 
fixed by the condition: a/ R{z) = z'^ + 0{z), z ^ oo. The branch assumes pure imaginary 
values of opposite sides on the edges of a, seen as a cut of C. 
Taking into account the equalities 

p{fi)dfi = 0, 



and 



0, / , . . = -TTZ, 



we find that 

p(A) = + ^v.p. / ^^Mv^^/^ . (A.15) 

The constant C can be found as follows. Denote f{z) the Stieltjes transform of p: 

Recalling that is real analytic, using p.8|l . and arguing as above, we find that 

fi,) = _ / - V'{z) dp 

2 27rz 7, p-z ' 

(the analog of the formula with Rq instead of R2 is valid for any finite number q of intervals 
in (in|)). 

Write now the minimum condition ()2.4j) as 

Veff{b)-V,ff{a) = 0, 

or, in view of ()2.6|) and ()A.16|) . as 

fi\) + ^^]d\ = 0. (A.18) 



The condition is valid for any potential, in particular, for V + eip. According to ()A.17j) 
for any sufficiently small e the integrand here is proportional to a/ -R2(A), in which a and b 
are now functions of e. Hence, the integrand vanishes at the edges of the support and the 
derivative of ()A.18|) with V replaced hj V + eip with respect to e at e = is 

/(A) + ^)rfA = 0. (A.19) 
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This and the formula 

■ ^ r KMi ^ 1 r m^/RM^^^ (^.20) 

following from ()A.15|1 . yield in the case v^(A) = tX: 

C P2(A) 
v/(62-A2)(A2-a2) 27rv/(&2-A2)(A2-a2) 

where P2(A) = A2 - {a'^ + h'^)/2, 

C=^, (A.21) 

' J -a v/(fe'-A2)(a2-A2)' ' ^(62_A2)(a2_A2)' 
Now we can find for v9(A) = tA. We have by fCTT]) and (fXT22|l : 

/3i := / p(A)dA = ;rV(^2Ji - /i J2), (A.23) 

where 

dA ^ P2(A)rfA 



By using standard formulas (see e.g. formulas (3.159)), we find 

/I - P \ 
h = 2K{k)/a, /a = 2a K{k) - E{k) , (A.24) 



Ji = ir(A;')/a, J2 = a (^P(A;') - l±^K{k')^ , (A.25) 

where K{k) and -E(/i;) are the complete elliptic integrals of the first and second kind, k = 
b/a, k' = Vl - A;2. These formulas and the identity EK' + E'K - KK' = 7r/2, where 
K' = K{k'), E' = E{k') (|2]], formula (8.122)) imply 

A more involved version of the above argument leads to ()3.22j) . We note first that to 
prove the formula for a real analytic ip it suffices to consider 

<^zoW = Y- — ' ^0 ^ ^■ 

A — Zq 

We have in this case (see below): 

where I is defined in ()3.23|) . 

Assuming that this formula is valid and using the Cauchy theorem to write 

1 f ^{zo)dzo 
27^^ Jc^ A - zo 
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where C^- is the contour encirchng a in the clockwise direction, we obtain in view of the 
hnearity of [3i[ip\ in ip: 



2m 



21 2m 



dzo- 



Co 



Now the relation y^R^(X^^ = -^/R2{X + iO) yields (HT^ . 

To prove ()A.26|) we use the general formulas ()A.15|) . ()A.20|) . and ()A.19|) with 



(A-zo)^ dzoX-zo' 

Arguing as in the case if{X) = tX above, we obtain an analog of ()A.23|) whose denominator 
contains I of (|3.23p instead of Ji of ()A.22|) . and whose numerator is a bilinear combination of 
integrals of \R2{n)\^^^'^ over [02, ^2] and of derivatives with respect to Zq of the integrals over 
[61,02] and [02,62] of [{f i-zo)\ R2{fi)\^^^]~^ and P2{i2)[{fi- zo)\R2{fi)\^^'^]~^ , where P2 is defined 
now by the relation ^/R^ = P2{z) + 0{1/z), z ^ 00. These integrals can be expressed 
via the complete elliptic integrals of the first, second, and third kinds. Furthermore, the 
complete elliptic integrals of the third kind can be expressed via the incomplete elliptic 
integrals of the first and the second kinds, whose arguments depend on zq. This allows us to 

obtain a formula for /3i[(^2q], whose numerator is expressed via the complete elliptic integrals 
of the first and the second kind and derivatives with respect to zq of the incomplete elliptic 
integrals of the first and the second kind. The derivatives are proportional to {R2{zo))~^^'^ 
(see [20., formulas (8.123)) This and a bit tedious algebra lead eventually to ()A.26|) . 

Another derivation of ()A.26|) is given in |H] (see formula (3.14) of the paper). The 
derivation is based on a two step procedure of minimization of ()2.3|) : the first step is the 
minimization over all unit measures with a given charge /3i e (0, 1) of the "band" [02, 62] of 
the support, and the second is the minimization of this minimum over 
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